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Abstract 

We introduce a novel approach to estimate the thermodynamic pressure from heavy-ion collisions 
based on recently measured higher-order moments of particle multiplicities by the STAR experiment. 
We start with htting the experimental results in the most-central collisions. Then, we integrate them 
back to lower ones. For example, we find that the first-order moment, the mean multiplicity, is exactly 
reproduced from the integral of variance, the second-order moment. Therefore, the zero-order moment, 
the thermodynamic pressure, can be estimated from the integral of the mean multiplicity, the possible 
comparison between such a kind of pressure (deduced from the integral of particle multiplicity) and the 
lattice pressure and the relating of Bjorken energy density to the lattice energy density are depending 
on lattice QCD at hnite baryon chemical potential and hrst-principle estimation of the formation time 
of the quark-gluon plasma (QGP). 

PACS numbers: 05.70.Ce, 05.45.-a, 25.75.Nq 
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I. INTRODUCTION 


The structure of matter at pressure and/or energy density much larger than the critical 
values defining the hadron-quark deconfinement phase transition(s), is to a large extend well 
understood. One of the yet-unsettled problems of theoretical physics is the characterization of 
the equation(s) of state (EoS) describing the behaviour of thermodynamic quantities at finite 
temperatures and densities. To this end, we utilize effective models, lattice QCD simulations 
and high-energy collisions. To the latter, the ultimate experimental evidence is coupled. On one 
hand, the lattice QCD calculations turn to be very reliable, especially at high temperatures and 
densities. On the other hand, the exact EoS of the hadronic matter is still rather complicated to 
be modelled. In describing the ground state properties of nuclear matter having large number 
of hnite nuc 


successful 


Q 


ei, the Hartree-Fock theories using Skyrme effective interactions are found quite 
3|. Nevertheless, it seems that serious concerns about basic physical symmetries 
arise, when using EoS derived from the Skyrme interactions j4|, especially at hnite temperatures. 

The early Universe is assumed to undergo a rapid phase transition from a phase dominated by 
colored deCTees of freedom (partons) to a phase dominated by color-neutral degrees of freedom 
(hadrons) |^. Such a phase transition is now routinely reproduced in heavy-ion collisions at the 
Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). The total number 
of produced particles exceeds several thousands. Therefore, one can expect that the generated 
system contains collective phenomenon. In light of this, we recall that the hydrodynamic 
description of the heavy-ion collisions goes back to Landau in early hfties of the last century. 
The hydrodynamics predicts long-range correlations induced by thermal noise. 

The statistical-thermodynamic quantities, like pressure, energy density and number density, 
are well-known tools to describe nature, degrees of freedom, decomposition, size and even 


fl. 


the overall dynamics controlling evolution of the medium from which they are originating 
The first-principle lattice calculations offer an excellent control on these quantities, especially 
at equilibrium. The effective models based on various approximations are also able to shed 
light on these quantities. The present work suggests a novel approach towards deducing the 
thermodynamic pressure from the high-energy collisions with different centralities. We would 
like to highlight that this estimate is approximative. As stated in the title, the proposed method 
is an estimate due the various sources of uncertainties, which will be discusses in sec. Ill Cl 
The success in estimating the thermodynamic pressure opens horizon for other thermodynamic 
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quantities, such as the energy density which can be compared with the Bjorken energy density. 
Based on the higher-order moments of particle multiplicity measured by STAR at RHIC, section 
m we propose to estimate the pressure. 

It is in order now to highlight possible procedure in estimating energy density, for instance. 
This is not a straightforward implementation of the proposed method. Estimating the ther¬ 
modynamic pressure which is directly related to the free energy, the energy density can be 
deduced from the derivative of the free energy with respect to inverse temperature. This obvi¬ 
ously illustrates the need to implement another variable different than the one responsible for 
the derivatives of the high-order moments, the chemical potential. In a future work will analyse 
the impact of this new variable. 

Reliable estimation of EoS (statistical-thermodynamic quantities) has deep impacts on var¬ 
ious cosmological and astrophysical aspects. For instance, viscous EoS deduced from lattice 
QCD and heavy-ion collisions was implemented in early Universe (QCD era) js 12|. Section 
m is devoted to the discussion and outlook. Our remarks about the proposed methis are 
elaborated in section III Cl 


II. HIGHER-ORDER MOMENTS OF NET-PROTON MULTIPLICITY 


The higher-order moments can be studied in different physical quantities. For example, 
the higher-order moments of charged-particle multiplicity distribution have been predicted four 
decades ago [l^. Recently, the higher-order moments have been reported by STAR at various 
RHIC energies and by lattice QCD calculations 19-21|. The empirical relevance of 


the higher-order moments to the experimental measurements has been suggested in Ref. 22j. 
Accordingly, the measurement of the correlation length ^ is very crucial. It is bounded by the 
hnite volume and lifetime of the system through causality principle. In a relativistic system, 
the growth rate of correlated domains should be ultimately bounded by c. 

The normalization of higher-order moments gives additional insights about their properties. 
From statistical point of view, the normalization is done with respect to the standard deviation 
a, which is related to i.e. r—th order moment is to be normalized to a^. Therefore, it 
provides with a tool relating moments with various orders to the experimental measurement. 
For instance, the susceptibility of the distribution, second order moment, gives a direct measure 
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for a. 

There are several techniques to scale the correlation functions. The survey system’s optional 
statistics module represents the most common technique, i.e. Pearson or product moment 
correlation. This module includes the so-called partial correlation which seems to be useful 
when the relationship between two variables should be highlighted, while effect of one or two 
other variables can be removed. The products of higher-order moments of conserved quantities 
can be directly connected to the corresponding susceptibilities and be directly related to long 
range correlations 
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25|. 


Some details about the experimental procedure are now in order. The STAR detector excel¬ 
lently identihes produced particles and has a large acceptance for the event-by-event fluctuation 
analysis. For each beam energy, an extensive quality assurance is performed in order to mini¬ 
mize the fluctuation of detector efficiency. From the Time Projection Chamber (TPC) center, 
the event-by-event net-proton distribution lying within ±30 cm along the z position of the 
interaction point and 2.0 cm radius in the transverse plane can be excellently detected. 

The geometry of each heavy-ion collision (and thus centrality of the collision) can be charac¬ 
terized by the number of participating nucleons (Npart)- By using centrality bin width correction 
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27|, hnite centrality bin width effect has been avoided. The charged produced particles 


can measured between 0.2 < pt < 2.0 GeV and \ri\ < 0.5, where pt and p are transverse mo¬ 
mentum and pseudo-rapidity, respectively. Centrality selection has been done by uncorrected 
charged particles measured within the wider window 0.5 < \p\ < 1.0 in order to avoid the 
auto-correlation effect in the higher moments analysis. 

The positive and negative charged particles distributions are assumed to be independently 
Poisson distribution where no dynamical correlation among the positive and negative charge 
particles are taken into account. Accordingly, the net-proton distribution is taken as a Skellam 


distribution 


28| . which is a discrete probability distribution of the difference between two 


statistically independent variables. Each of these random variable has Poisson distri 


a different expected value. The Skellam distribution is baseline for this analysis 29 


Dution and 
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Fig. 1 ; Mean (left panel) and standard deviation (right panel) are plotted as function of the baryon 


chemical potential The experimental results (symbols) are taken Ref. 
Eqs. ([2]) and (l3|), respectively. 


29l |. The curves represent 
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Fig. 2 : Skewness (left panel) and kurtosis (right panel) are plotted with 


respect to the baryon chemical 


potential The experimental results (symbols) are taken from Ref. 29(]. The curves represent Eqs. 
([H) and ([5|), respectively. 

A. Integrals of higher-order moments 

The experimental higher-order moments at the largest (Npart) are related to most-central 
high-energy collisions and likely comparable to lattice QCD simulations. The effects of varying 
centralities of the collisions shall be discussed in section III Dl The first four moments of net- 
proton multiplicity are plotted with the baryon chemical potential fib, which is related to the 
nucleon-nucleon center-of-mass energy . For simplicity, we write y/s instead of ^/SNN■ 

1.31 ±0.03 [GeV] 

“ 1 ± (0.27 ± 0.01 [GeV^D^sVV ^ ^ 

We notice that M and S increase with increasing fib, while the standard deviation a decreases 
and kurtosis k remains almost constant. The given curves represent the statistical fits, Eqs. 
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(12])-([5]). The baryon chemical dependence of mean, standard deviation, skewness and kurtosis, 
respectively, reads 


M — 

+ C2, 

S = + C 3 , 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


K = 0.004 ±0.002 


where Oi = 37.8±0.11, h = 0.6±0.11, Ci = -2.105±1.1, 02 = 175.75±11.07, h = -5.27±1.02, 
Cl = 89.5 ± 13.13. as = 0.16 ± 0.06, 63 = 1-5 ± 0.5 and C 3 = -0.001 ± 0.0005. Equations dH)- 
(15|) are analytical expressions representing the first four orders of moments. The procedure 
introduced in the present work is the integration of the moments of particle multiplicity. This 
consecutively results in moments with lower orders. By implementing this we propose that the 
thermodynamic pressure can be estimated. 

B. Thermodynamic pressure in high-energy collisions 

Our approach towards deducing thermodynamic pressure from the higher-order moments 
measured in the heavy-ion collisions is based on: 

• the higher-order moments of particle multiplicities are accessible in the high-energy ex¬ 
periments through the correlation length 

• fitting the experimental results on higher-order moments makes it feasible to integrate 
them back to lower ones, and 

• therefore, the zero-order moment of mean particle multiplicities, the intensive thermody¬ 
namic pressure, can be obtained, accordingly. 

In order to test the given approach, we can apply it on the variance distribution, for example. 
The integral of the variance gives an estimation for the mean multiplicity. In Fig. [31 the STAR 
results on normalized mean multiplicity (symbols) are given as function of the baryon chemical 
potential, fXb. The translation between y/s and /i^ is governed by Eq. ([1]). The integral of 
variance, Eq. (|3|) and left-hand panel of Fig. [H is represented by the thick solid curve. 



( 6 ) 


7i 
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Fig. 3 : The STAR results on normalized mean multiplicity (STAR BES I) are given in dependence 
on the baryon chemical potential, The cnrve represents the integral of Eq. ([3|), left-hand panel of 

Fig. m 


where ai = 25.94 ± 3.8, /3i = 58.6 ± 3.7, 71 = —5.27 ± 1.024 and = 11. It is quite obvious 
that the two curves in Fig. [3] and left-hand panel of Fig. [1] are likely identical. This makes 
it eligible to integrate the mean multiplicity, Eq. ([2]) right-hand panel of Fig. [U in order to 
deduce the thermodynamic pressure. Fig. 01 As expected, the overall behavior is monotonic. 
Increasing leads to an increase in the resulted pressure. This behavior was observed in lattice 
QCD at hnite fib [Mt]. When fib is translated into the dependence of integral of mean or 
double integrals of variance (dimensionless pressure), Eq. ([6]), on -y/s is given in right panel. 
Here, the behavior is non-monotonic. It is worthwhile to notice that the relation fib{\fs) was 
phenomenologically estimated in the hnal state of particle production. 

At small yi, the resulting pressure is large. This is obviously related to the large stopping 
power related to the incident energy. When the value of ^/s reaches about 1 — 2 GeV, the 
pressure considerably decreases. This would be related to a drastic change in the effective 
degrees of freedom. At large the resulting pressure gets relatively small values. Further 
raise in ^/s does not change the thermodynamic pressure. The overall behavior can be described 
as ’z-shape’ or nearly ’tangen-hyperbolicus’-like. 

The integral of mean multiplicity (dimensionless thermodynamic pressure) reads 


p{fib) = a^fib H-+ ^5, 

C5 + i 

T(V^) = exp (c6\/s) + d^y/s, 

I Cq 


(7) 

( 8 ) 
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where 05 = -2.1±1.2, h = 37.8±1.97, C 5 = 0.6±0.011, rfs = H, ae = 89.5±1.3, h = 58.6±3.9, 
Ce = —5.27 ± 0.02, and = 11. This result is obviously original. So-far this quantity was not 
accessible in the high-energy experiments. An extension to other thermodynamic quantities is 
therefore planned in future works. 




sqrt(s) [GeV] 


Fig. 4: The integrals of mean multiplicity distribution is given in dependence on the baryon chemical 
potential (left) and center-of-mass energy (right). The vertical bands show the critical deconfinement 
phase transition. The overall behavior has is nearly ’tangen-hyperbolicus’-like. 


To make the new method more clear we summarize the procedure, shortly: 

• From consecutive integrals of the high-order moments, an analytical expression for the 
number density (the multiplicity) can be deduced. 

• By integrating this expression with respect to the chemical potential, the thermodynamic 
pressure will be estimated. 


C. Remarks on proposed estimate 

We want to highlight that the present work does not claim that this is an exact method to 
deduce the thermodynamic pressure. As we discussed earlier, to this end, we still have to take 
into consideration various uncertainties. 

• First, the direct comparison of the pressure (deduced from integral of mean multiplicity, 
especially in central collisions) with the lattice pressure is not possible. 

• Second, the comparison of the Bjorken energy density with the lattice energy density is 
also not allowed. An estimate of the thermodynamic energy density from the high-energy 
collisions seems to implement sophisticated methods. 
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• Third, we should verify the assumption that the moments of net-protons are almost 
equivalent to the ones of all baryons should be correct. Nevertheless, the assumption 
that the net-proton higher-order moments are nearly equivalent to that of net-baryon 
finds supporters [32 1. 

• Fourth, seeking for completeness, we recall the electric charge distribution. It is conjec¬ 
tured that their moments are related to the Taylor expansion coefficients of the pressure 
with respect to electric charge chemical potential at the freeze-out. This would be insuf¬ 
ficient to determine the pressure, since the overall integration constant cannot be fixed. 


• Fifth, the experimental measure for the higher moments utilizes particle multiplicity and 
very rarely the correlation length. Thus, the possible highlights on the critical phenomena 
related to hadron-quark phase transition(s) are apparently limited to statistical insights. 


In relating our estimate of the thermodynamic pressure to that deduced from lattice QCD 
simulations, the difficulty of performing lattice calculations at finite baryon chemical potential 
fib should be taken into account. Assuming chemical freeze-out |45l-l48| . the latter is to be 
related to the center-of-mass energy ^/s. 

In relating our estimate of the thermodynamic pressure to the hydrodynamics, the formation 
time r should be measured. The guess that r varies with ^/s should be first verified, before 
confronting and/or merging pressure and energy density measured in high-energy collisions and 
calculated in lattice QCD. So far, there is no direct measurement or first-principle estimation 
for the QGP formation time r. In a very simple treatment, the formation time of partons is 
indirectly related to their average energy The hadrons are emitted only from the surface 


of freeze-out 


45 


48l |. whereas direct photons are conjectured to be emitted throughout the life¬ 


time of the evolving system. Therefore, electromagnetically and/or weakly interacting probes 
would carry information about the formation time. For the QGP temperature, we have so-far 
two different measurements using direct photons. Further measurements are very desired. Also 
the estimation of the unified energy density and its relation to the critical temperature using 
traverse mass of different particles should be extended. 
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D. An estimate of the higher moments 


Fig. |5] shows the dependence of the hrst four moments on centrality of the collision, which 
is geometrically dehned by the averaged number of interacting nucleons (Npart)- As given in 
previous sections, the maximum (Npart) is likely related to the most-central collisions. Varying 
(Npart) from 50 to 200, we move from peripheral to central collisions. We notice that increasing 
centrality (here marked by vertical arrow) is accompanied by increasing mean and variance and 
by decreasing skewness and kurtosis. The solid curves represent the hts of the most-central 
results as given in Figs[T]and[2l Only for mean and variance, the hts for different centralities are 
shown. From these two quantities, we want to check again the validity of the given approach 
and then deduce the variation of the thermodynamic pressure with the collision centrality. 





tib [GeV] 


[GeV] 


Fig. 5 : Mean (a), variance (b), skewness (c) and kurtosis (d) are given dependence on the collision 
centrality. Increasing centrality (marked by vertical arrow) is accompanied by increasing mean and 
variance and by decreasing skewness and kurtosis. Solid curves present the fits of the most-central 
results as given in Figs [T] and [2j For mean, the fits for different centralities are shown. 


Left-hand panel of Fig. |6] shows integral of mean at various centralities in dependence on 
the baryon chemical potential. As expected, the overall behavior is monotonic. Reducing 
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Fig. 6: Left-hand panel: integral of mean multiplicity at various centralities is given in dependence on 
the baryon chemical potential. Right-hand panel: double integrals of variance at different centralities 
are given in dependence on center-of-mass energy. Assuming that the critical deconfinement phase 
transition does not depend on the centrality, one vertical band referring to that is drawn. 


the centrality reduces the values of integrals of mean (pressure). Right-hand panel shows the 
dependence of double integrals of variance at various centralities in dependence on center-of- 
mass energy. Here, the behavior is non-monotonic. The overall behavior can be described as 
follows. At small i/s, the resulting pressure is large. The large stopping power is likely the 
reason. When ^/s reaches 1 — 2 GeV, the resulting pressure considerably decreases. Assuming 
that the critical deconhnement phase transition does not depend on the centrality, one vertical 
band referring to that is drawn. A drastic change in the effective degrees of freedom is likely. 
At large ^/s, the resulting pressure gets relatively small values. Further raise in y/s does not 
change the pressure. 


III. DISCUSSION 

Current particle detection technologies are mainly based on properties of the leptons and 
hadrons, for instance. In high-energy experiments, the inputs are hadrons (or leptons) acceler¬ 
ated to high energies and the outputs are produced colorless particles which hnally freeze out, 
chemically and thermally, and then will be detected. The running strong coupling prevents 
quarks to escape from conhned hadrons. These are the experimental obstacles against direct 
or even indirect detection of thermodynamic quantities in early stages of high-energy collisions. 
Fig El 

The status of estimating the physical properties of hadronic and partonic matter is very 
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Fig. 7; Schematic illustration of various stages of a central heavy-ion collision. Measured and calcu¬ 
lated thermodynamic quantities in partonic and hadronic phases are outlined. 

shortly summarized in Fig. [71 The first-principle lattice QCD calculations shed light on ther¬ 
modynamics of the strongly interacting matter. Not only on numerical techniques the lattice 
calculations are depending, but also on the assumption of equilibrium statistical thermody¬ 
namics. The latter refers to a system that evolves to a state which is independent of the past 
history and time. The equilibrium state is specihed by energy, volume and number of particles. 
For a system at equilibrium, there exists a positive differential entropy function which increases 
with energy at hxed volume and particle number. In other words, the entropy is only dehned in 
an equilibrium system. Non-equilibrium phenomena, such as heat conduction or diffusion, are 
apparently of great interest in high-energy physics. The non-equilibrium statistical mechanics 
distinguishes fast variables (characterized by microscopic time and space scales) from slow ones 
(characterized by macroscopic time and space scales) [34 1. 

The physics program of high-energy experiments is very successful in deducing essential 
physical quantities like identihcation of produced hadrons, measuring their physical properties, 
counting their multiplicities, and estimating their fluctuations and correlations in different 
rapidity intervals, etc. The hadron multiplicities, fluctuations and correlations are experimental 
tools to analyse nature, composition, and size of medium from which they are originating. 
Based on analysing the particle abundances or momentum spectra, the degree of equilibrium 
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of the produced particles can be estimated. The particle abundances can help to establish the 
chemical composition of the system. The momentum spectra can give additional information 
on the dynamical evolution and the collective flow. For instance, at chemical equilibrium, the 


entropy and number of produced particles are hxed. Thus, statistica. 
be applied in order to estimate various thermodynamic quantities 


the lattice calculations was found excel 
Therefore, the QCD phase diagram 


ent 


40 


35 


^^ermal approaches can 
. The agreement with 


44| . Both are likely related to central collisions. 


44l | and freeze-out curves can be mapped out |45l448| 


IV. SUMMARY AND OUTLOOK 

Based on the higher-order moments of particle multiplicities, which are experimentally ac¬ 
cessible through the correlation length We started with htting the experimental results, we 
introduced an approach to estimate the thermodynamic pressure from the heavy-ion collisions. 
Afterwards, we integrated them with respect to the chemical potential. This resulted in the 
lower-order moments. For example, the hrst-order moment, which is the mean multiplicity, 
was exactly reproduced from the integral of the variance. Fig. [31 which is the second-order mo¬ 
ment. Accordingly, the zero-order moment, which is noting but the thermodynamic pressure, 
was obtained from the integral of mean multiplicity. Fig. |4l The concerns about the certainties 
of the proposed method have been elaborated. 

One of the ultimate goals of the comprehensive program of high-energy experiments is a 
reliable estimation for all equation of state and that of motion. So far, we measure collective 
properties of produced particles, such as multiplicities, correlations and fluctuations. Further¬ 
more, the trustful order-parameters characterizing the phase transition are related to thermo¬ 
dynamic quantities. The lattice QCD simulations at hnite temperature and density make use of 
such order parameter in determining the critical phenomena accompanying the deconhnement 
to the colored phase, the quark-gluon plasma. In the present, we introduce an estimate of 
the thermodynamic pressure from the measured multiplicities. Should the proposed results be 
complete and accurate, we shall be able to compare the experiments directly with the lattice 
QCD simulations. Obviously, this brings the high-energy experiments to a very bright horizon 
and enables tackling various still-unsolved problems, for example, orders of the deconhnement 
phase-transition and critical temperature and exponents. The latter are related to the uni- 
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versality classes of the system of interest. For an accurate estimation of the thermodynamic 
quantities: 

• quite gauge of the phase space, 

• comprehensive measure of the baryon multiplicity, 

• full identihcation of the produced particles, 

• exact determination of the hreball volume, 

• complete estimation of possible correlations and fluctuations, 

• characterizing the geometrical (spacial) and temporal evolution of the hreball, 

• complemented metering of the freeze-out boundaries 

• etc. 

are crucially signihcant. No doubt that advanced experimental facilities are indispensably 
essential to fulhl these demands. 
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